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ABSTRACT 

Magnetic Bright Points (MBPs) are among the smallest observable objects on the solar photosphere. 
A combination of G-band observations and numerical simulations is used to determine their area 
distribution. An automatic detection algorithm, employing 1-dimensional intensity profiling, is utilized 
to identify these structures in the observed and simulated datasets. Both distributions peak at an 
area of «45000 km 2 , with a sharp decrease towards smaller areas. The distributions conform with log- 
normal statistics, which suggests that flux fragmentation dominates over flux convergence. Radiative 
magneto-convection simulations indicate an independence in the MBP area distribution for differing 
magnetic flux densities. The most commonly occurring bright point size corresponds to the typical 
width of intergranular lanes. 

Subject headings: Sun: evolution — Sun: granulation — Sun: photosphere — Sun: surface magnetism 



1. INTRODUCTION 

The dominant pattern covering the quiet solar pho- 
tosphere is the granulation. Plasma flows remove the 
magnetic flux from granules into the dark inter-granular 
lanes where it clusters together to form small magnetic 
concentrations of 1" or less, with field strengths often in 
excess of a kilogauss. Magnetic Bright Points (MBPs) 
are a manifestation of these kilogauss fields. They are 
among the smallest observable objects on the photo- 
sphere appearing as intensity enhancements within the 
inter-granular lanes. They are formed by complex pro- 
cesses involving the interaction of magnetic fields with 
the convectively unstable plasma and can provide a con- 
duit for channeling kinetic en ergy into the upper atmo- 
sphere (|de Wiin et al. 1 120090 . Although all MBPs are 
found to reside in areas of high magnetic field, the oppo- 
site is not always true, with some high m agnetic field in 
plage regions appearing devoid of MBPs (Ishikawa et al. 
120071 ). An examination of the physical parameters re- 
quired for a MBP to exist, i.e. minimum or maximum 
size, and lower magnetic field threshold, will further our 
knowledge on the creation and evolution of such small- 
scale photospheric magnetic fields. 

The physical processes re sponsible for MBP fo rma- 
tion have bee n simu lated by iSchussler et al.1 ()2003l ) and 
IShelvag et al.l (|2004f ) . Magneto-convection models for 
convectively unstable plasma in the photosphere, are 
combined with the radiative properties of that plasma 
and include the partial ionization of hydrogen and the 
other most abundant elements. The validation of numeri- 
cal simulations by observations, may allow us to conclude 
whether the physics in the simulations can describe the 
real Sun and hence help us interpret the observational 
findings. 

G-ban d imaging forms a common bas is for MBP obser- 
vations. ISanchez Almeida et al.l (|2004f ) find that MBPs 
have a maximum intensity 1.8 times the mean pho- 
tospheric value. Several authors propose that such 
intensity enhancements are caused by a significant 



weakening of the CH absorption lines, due to the 
dissoci ation of CH molecules at i ncreasing tempera- 
tures (ISanchez Almeida et al.l [200ll iSteiner et all 120011 : 
lLanghans et al.l I2002D . The latter group identify two 
types of bright points that occur in G-band images; those 
that are associated with magnetic structures, and others 
that exist at the edge of rapidly expanding granules. The 
latter are believed to be density enhancements, caused by 
a build-up of material attempting to flow down the thin 
inter-granular lanes at the perimeter of a granule. 

The sheer number of MBPs requires automated al- 
gorithms for their detection and tracking. In general, 
higher s patial resoluti o n lead s to smaller sizes being de- 
tected. iBerger et al.l (|1995f ) applied an altered blob- 
finding algorithm to separate MBPs from granules. They 
quantified the MBP size in terms of the FWHM inten- 
sity diameter, taking the smallest dimension across the 
identified objects. A semi-automatic process was imple- 
mented that included nonlinear least-squares Gaussian 
fitting to the observed intensity profiles. Their analysis 
revealed a modal diameter of 220 km, an average diam- 
eter of 250 km, a nd a diameter range of 120-600 km. 
iBovelet fc Wiehrl ([2003D employed a specialized version 
of the multiple-level tracking pattern recognition soft- 
ware (jBovelet fc Wiehrl 120011 ) . which applies several de- 
creasing intensity thresholds to an image. A dominant 
diame ter of 220 ± 25 km was detected. iWiehr et al.l 
(2004) repeated this procedure using higher spatial res- 
olution observation s, and found a predominant dia meter 
of 160 ± 20 km. ISanchez Almeida et all (|2004D visu- 
ally identified MBPs in individual G-band images. The 
area was determined by a segmentation algorithm, and 
their diameter measured by fitting a double gaussian pro- 
file. They obtain 135km for the minor axis, which seems 
to be set by the angular resolution of the observations. 
This double gaussian decontaminates the profiles, by tak- 
ing into account the intensity distribution of the dark 
local background, wit hi n which the MB Ps are situated 
(|Title fe Bergedfl99ll . lUtz et al.1 (pool utilized an al- 



2 



P. J. Crockett et al. 



tered version of the iBovelet fc Wiehrl ([2001D algorithm, 
and applied it to Hinode SOT observations. The size of 
MBPs was defined by placing an upper and lower inten- 
sity threshold on the segmented structures, resulting in 
mean diameters of 166-218 km. 

In this paper we use observations and numerical simu- 
lations to investigate the area distribution of MBPs. An 
au tomatic detecti o n and tracking algorithm, described 
in iCrockett et al.l f|2009T ) (hereafter, Paper 1), is devel- 
oped further and applied to high resolution G-band im- 
ages. Section 2 discusses the observations, with emphasis 
on an automated algorithm used for MBP detection and 
size determination. A description of the numerical sim- 
ulations are given in Section 3. Our main findings are 
presented in Section 4, with concluding remarks in Sec- 
tion 5. 

2. OBSERVATIONS AND DATA ANALYSIS 

The data was obtained by the newly-commissioned 
Rapid Oscillations in the Solar Atmosphere (ROSA) in- 
strument, installed at the 76 cm Dunn Solar Telescope 
(DST), in New Mexico, USA (IJess et al.ll20loft . The ob- 
servations were taken on 28 May 2009 through a 12 A 
filter, centered at 4305A (G-band), during a period of 
excellent seeing. Post-facto s peckle reconstruction al- 
gorithms (W oger et al. I [20081) . in addition to rigorous 
image de-stretching using a 40 x 40 grid, (e quating to 
a sal .7" separation between spatial samples, ([Jess et al.l 
2008)), was implemented to remove effects caused by at- 
mospheric seeing. We observed a 70" x 70" quiet solar 
region at disk center, achieving diffraction-limited imag- 
ing with 0.069" pixel" 1 . Figure [T] displays a typical G- 
band image from the dataset, with multiple MBPs visi- 
ble in the central region. Analysis of the data was per- 
formed with an updated version of a detection algorithm 
described in Paper 1, which uses intensity thresholding 
to map the intergranular lanes. A compass search al- 
lows MBPs to be disentangled from bright pixels within 
granules, while object growing accounts for any pixels 
that might have been removed when mapping the lanes. 
One of the disadvantages of the algorithm described in 
Paper 1 is the requirement for the image to be divided 
into segments, with each subsequent segment being pro- 
cessed individually. Here we use an updated algorithm 
which operates on the entire 70" x 70" image sequence, 
thus improving computational time and accuracy. This 
development is particularly important, as it permits ac- 
curate estimates of MBP areas. Mapping the location of 
the inter-granular lanes, with an overestimation of the 
intensity threshold, is used to separate out bright ob- 
jects. The threshold set is the mean image intensity plus 
one sigma. All structures under this level are considered 
a lane and are not investigated by our algorithm. The 
vast majority of MBPs retain higher intensities, however 
some very dull MBPs may be lost at this stage. Objects 
are then investigated individually. We impose a 3-sigma 
intensity variation limit on each object, in order to fully 
separate MBPs from the granules. Any bright object 
possessing an intensity range greater than 3 sigma is bro- 
ken up into smaller objects, until the resulting structures 
comply with this condition. 

The detection of MBPs is carried out by an extended 
version of the compass search (see Paper 1 §4.3), and in- 
corporates gradient thresholding through intensity pro- 



filing. A one-dimensional variation in intensity, across a 
selected region of the image, is first determined (see left 
panel of Fig. [2]). Intensity profiles for each object are es- 
tablished for eight directions, symmetrically positioned 
about the objects centre-of-gravity. The stipulation that 
a lane must be in close proximity to the MBP remains. 
The algorithm now actively searches for inter-granular 
lanes, by using the turning points of the intensity pro- 
files, which arc located at the centre of the lanes (left 
panel of Fig. [2|). Hence, the lanes are located from in- 
situ intensity profiling. Each measurement is specific, 
not only to individual objects, but in every considered 
direction as well. To establish each turning point the 
one-dimensional line, from which intensity profiles are 
procured, is extended until two stationary points exist in 
the profile, i.e. where the rate of change in intensity (y) 
as a function of distance (x) is zero, dy/dx = 0. A limit 
on the distance between the turning points eliminates 
large objects, such as granules. 

Gradient thresholding is applied to all intensity pro- 
files to disentangle MBPs from granules. MBPs retain a 
very steep intensity change in all directions, compared to 
a more gradual variation associated with granules. The 
maximum gradient is determined from any part of a pro- 
file falling between the two turning points (see Fig. [5]) . 
The threshold gradient is derived for each individual im- 
age through a selection of 500 random objects, and is 
calculated prior to the execution of the algorithm. A 
threshold is determined by adding a 1-sigma value to the 
median gradient recorded for each image. 

A significant improvement of the present algorithm 
concerns the growing of MBPs. A newly developed pro- 
cess provides an independent threshold range for each 
object for accurate area representation. The algorithm 
rotates a one-dimensional line through 360 degrees, in 5 
degree steps, about an object's centre of gravity. Inten- 
sity values at the turning points of the profiles, i.e. the 
lanes, are acquired at each angle. To aid the accurate 
determination of turning points, the data is re-binned 
by a factor of ten and smoothed. Thus a narrow inter- 
granular lane and the associated turning point can be 
clearly identified. The maximum turning point intensity 
is taken to provide a lower cutoff to our growing algo- 
rithm, while the upper boundary is set as the maximum 
intensity level occurring within the seed region. The 
growing procedure includes any conjoining pixels that 
are above the lower threshold cutoff. The MBP area is 
determined by totaling the number of pixels within each 
structure. Our sampling of 0.069" pixel -1 provides an 
area of 2500 km 2 pixel -1 . This procedure, demonstrated 
in Figure [3[ reproduces 90% of MBPs to within a 10% 
error of visually identified areas. Setting the lower inten- 
sity threshold as the brightest surrounding lane enforces 
an upper limit on the area of the MBPs, i.e MBPs are 
grown to their maximum dimensions. 

3. NUMERICAL SIMULATIONS 

We use the MURaM code ([Vogler et al.ll2005h to carry 
out simulations of radiative magneto-convection in the 
upper solar convection zone and photosphere. This code 
uses a fourth-order, central difference scheme for com- 
puting the spatial derivatives, and a fourth-order, Runge- 
Kutta scheme to advance the solution forward in time. 
The solution is stabilized against numerical instabilities 
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using additi onal artificial hyperdiff u sive terms, de s cribed 
in d etail bvlCaunt fc Korpil (l200l. IVorier et al.l (|2005h 



and IShelvag et al.l (|2008l) . The size of the computational 
domain used for the simulations is 12 x 12 x 1.4 Mm 3 , 
resolved by 480 x 480 x 100 grid cells providing a res- 
olution of 25 km per grid cell. However, we emphasize 
that as a result of hyperdiffusivity, the size of the small- 
est structures produced in the smulations can be larger 
than a single grid cell. Due to the dependence of the hy- 
perdiffusivity coefficients on the local solution, it is not 
possible to globally define a quantity, uniquely represent- 
ing the resolution. However, a standard test, such a s the 
stron g (compression ratio 100) Riemann shock tube (jSodl 
can be used to provide an indication of the res- 
olution. The results of such tests for similar codes (i.e. 
4-6-th order central difference spatial scheme and hyper- 
diffusive sources) show that even for such an extreme 
case, the shock front is diffused over 2-4 grid cells, de- 
pending on the relative position of the shock front with 
respect to the grid (Caunt & Korpi, 2001; Shelyag et al., 
2008). Consequently, the resolution of the code for this 
case is about 50-100 km, a value similar to the resolution 
of the observations. The side boundaries are periodic, 
the upper boundary is closed for vertical and stress-free 
horizontal plasma motions, while the bottom boundary 
is transparent. The level corresponding to the visible so- 
lar surface is located approximately 400 km below the 
upper boundary. This setup allows us to perform radia- 
tive diagnostics of G-band images, and directly compare 
them with the observations. A detailed descrip tion of 
the method used is given in IShelvag et aP d200l . Here 
we provide a brief description of the process. 

For each of the light rays corresponding to a vertical 
plasma column in the simulation, we compute the LTE 
spectrum in the 4295-4315 A range, which consists of 
328 absorption lines, 239 of which are produced by CH 
molecules. The calculated spectrum is convolved with 
the G-band filter function. The magnetic splitting of CH 
lines, and its influence on G-band intensities, are suf- 
ficiently small for this effect to be neglected. G-band 
images obtained using this technique reproduce the dy- 
namic and radiative properties of magneto-convection, 
and show a large number of G-band bright points, cor- 
responding to the heated, and partially evacuated, mag- 
netic flux tubes seen in Figure [TJ 

4. RESULTS 

A series of 500 images were investigated, incorporating 
a total of 63312 MBPs. The MBPs cover approximately 
0.42% of the solar surface with a variance between 0.33% 
and 0.53% across the time series. Figure [4] displays the 
area distribution of MBPs, with their occurrence nor- 
malized to the mean number detected across all bins. 
The distribution was created by a summation of MBPs 
across all images, in 1 pixel bins. This technique may 
lead to the "double counting" of MBPs, some of which 
may have longer lifetimes than others. However, snap- 
shots of single frames produce a similar distribution, with 
approximately the same peak. Therefore, this technique 
is equivalent to the summation of multiple snapshot dis- 
tributions, each with similar parameters, resulting in an 
overall identical distribution. We therefore believe that 
any double counting does not pose a problem in our in- 
terpretation. 



The nature of the distribution appears to conform with 
log-normal statistics. To confirm this, a log-normal prob- 
ability density function (PDF) of the form, 
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where fi and a are, respectively, the mean and standard 
deviation of In x, is fitted to the data. Values of ^=3.25 
and cr=0.65 produce an excellent fit, shown by the over- 
plotted solid red line in Figure [4] To qu antify the good- 
ness of the fit, the x 2 error statistic (Wa ll fe Jenkins! 
2003) of the form, 
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where Oi and Ei are, respectively, the observed and ex- 
pected frequencies, is utilized. The observed frequencies 
correspond to the values obtained in the data, whilst the 
expected frequencies correspond to the theoretical val- 
ues set by the log- normal fit. Comparison of the real 
data with the fitted distribution reveals a conformity of 
99.5%, and confirms the MBP area distribution is well 
described by log-normal statistics. 

The peak of the distribution occurs at an area of 
45000 km 2 . Assuming a circular geometry, this corre- 
sponds to a diameter of 230 km. While this estimate ap- 
pears in general agreement with earlier works, there are 
a num ber of points that need to be emphasized lUtz et al.l 
([20091) find diameters of 218±48 km using a spatial sam- 
pling of 0.108" pixel -1 on Hinode SOT but this is depen- 
dent on the spatial sampling. Reducing the spatial sam- 
pling to 0.054" pixel -1 gives a diameter of 166±31 km. 
The latter value i s in agreement with the results of 
IWiehr et al.l ((2004) . Differences in the diameters may 
be explain ed by differences in the detection algorithms 
employed. IWiehr et al.l (|2004| ) employ the Multi Level 
Tracking (MLT) algorithm which utilizes decreasing in- 
tensity levels to identify and separate objects. MLT sets 
an initial uppermost intensity level. Bright structures 
which exceed this threshold are tagged. The intensity 
level is then lowered. Pixels above the new intensity 
level, adjacent to the structures identified in the pre- 
vious level, are added. New structures that appear at 
this level are tagged separately. This repetitive proce- 
dure is terminated after a last extension to a final in- 
tensity level deemed adequate for representing the ob- 
served pattern. The structures are forced to be separated 
by 2 pixels on all sides. The enforced separation and 
the somewhat arbitrary final threshold level, may affect 
the dimensio ns of the MBPs measured by missing dim 
edge pixels. lUtz et all (|2009f) employ a similar repeti- 
tive intensity thresholding technique to separate granules 
and MBPs. They impose an upper and lower intensity 
boundary on the pixels of an object to determine its size. 
The upper boundary is given by the maximum intensity 
in the object whilst the lower boundary is defined as the 
maximum minus 30% of the mean photospheric intensity. 
Again these conditions may limit the final size of MBPs 
to the brightest pixels. This effect may be exaggerated 
at the higher spatial sampling. Instead our algorithm 
sets a threshold which is specific to each MBP. Moreover 
by taking the threshold as the highest intensity within 
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the lane, we determine an upper limit in the MBP area 
including dim edge pixels which may be neglected in the 
earlier studies. 

In Figure IU the vertical dashed line at 10000 km 2 
marks a 4-pixel (2x2) resolution threshold, thus placing 
a limit on the smallest structures that can be resolved in 
our observations. The procedure of identifying the small- 
est structures in the observations was tested by convolv- 
ing the simulated G-band images with the Airy function 
corresponding to a 76 cm aperture and rebinning them 
to the spatial sampling of the observations. This test 
showed that 11 out of 12 MBPs, each covering an area of 
3-5 pixels in the degraded images, correspond to MBPs 
in the original (non-degraded) images. The sharp drop 
in the area of MBPs below the peak, implies that un- 
derlying physical processes are limiting the creation, and 
evolution, of very small structures. In addition, the dis- 
tribution shows only a relatively small number of large- 
scale MBPs with area greater than 200000 km 2 . 

Errors in the distribution were determined by compar- 
ing the visual estimates of the area with the output of 
the algorithm. The comparison reveals the distribution 
of errors for each bin from which a sigma value was deter- 
mined. At the peak of our distribution the error is ±8% 
(±3600 km). The uncertainty at the high end of the 
distribution decreases whilst, as we approach the diffrac- 
tion limit the errors become larger due to small number 
statistics. 

The simulations allow us to study the effect of chang- 
ing the net magnetic flux density on solar granulation 
and how this modifies the properties of MBPs. To pro- 
duce a theoretical distribution of their area, we first com- 
puted a series of consecutive G-band snapshots, based 
upon different average magnetic flux densities. The same 
procedures for MBP detection and area estimation were 
then applied directly to these simulations. We utilized 
100, 200, and 300 G signed vertical magnetic flux simula- 
tions, and compared the resulting distributions directly 
with the observations. The spatial resolution of simu- 
lations has not been degraded (Fig |4}. All three dis- 
tributions, based on different initial magnetic flux den- 
sities, agree with the log-normal form of the observed 
MBP area distribution. Crucially, these simulated G- 
band images also exhibit the same sharp decline from 
the peak to the diffraction limit. The peak of the distri- 
butions do not alter significantly, with a simulated maxi- 
mum at ~50000 km 2 . The simulations were degraded by 
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re-binning the data to match the observational scale of 
50km per pixel. A further comparison revealed that the 
simulated distribution remained unchanged. 

5. CONCLUDING REMARKS 

We utilize an improved automatic detection algorithm 
to study the area distribution of solar magnetic bright 
points (MBPs). The area distribution of MBPs follows 
log-normal stat istics. We interpret this result in a sim- 
ilar fashion as iBogdan et al.l ([1988), who suggest that 
the underlying fragmentation process is responsible for 
log-normal distributions. Similarly, we can not rule out 
coalescence. 

The peak of our MBP distribution occurs at 45000 km 2 
significantly higher than our telescope diffraction limit, 
and is consistent with the results of the radiative MHD 
simulations. The minimal area of MBPs is most likely 
defined by the width of the intergranular lanes, which 
is subsequently limited by the radiative and convec- 
tive energy balance, and mass conservation in magneto- 
convective processes. The peak in the area distribution 
seems to correspond to the most probable width of the 
intergranular lane. The area of large MBPs may be lim- 
ited by the lack of sufficient radiative heating in the 
larger flux tubes. As has been demonstrated (see e.g. 
iBerger et al.lll995T) , wall heating of magnetic flux concen- 
trations is not sufficient to increase the vertical radiative 
flux for flux tubes greater than 500 km in width. Thus, 
MBPs cannot be generated in large diameter magnetic 
flux tubes. Elongated bright points of large size may 
still, in principle, be formed. However, their formation 
will be inhibited by strong plasma motions and granule 
fragmentation. 
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Fig. 1. — Left: A 70" X 70" G-band image from the ROSA dataset. Middle: A 12 X 12 Mm section of the G-band image. Right: G-band 
simulation for an average field of 200G. 
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Fig. 2. — Top: Intensity profile typical of a MBP. A steep intensity gradient is followed by two minimum turning points at the centre 
of the dark lane. The algorithm detects the minimum turning points, marked by dashed red lines, in 8 separate directions. Bottom: The 
white marks indicate the locations of turning points in the image. These turning points reside in the center of the inter-granular lanes. 
The MBP shown has an area of 127,500 km 2 . 
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Fig. 3. — A line is rotated, in 5 degree steps, around an object to create one-dimensional intensity profiles (right). The turning points 
are marked by dashed red lines, whilst the green and blue dot-dash lines indicate the intensity levels of the left and right turning points, 
respectively. The maximum turning point, across all profiles, is imposed as the lower threshold for MBP growing. Note the re-binning of 
the data to ensure turning points are located accurately. 
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Fig. 4. — Top: The observed area distribution of MPBs. The green and blue dashed lines mark, respectively, the peak of the distribution 
at 45000 km 2 , and the diffraction limit at 10000 km 2 . A log- normal fit to the distribution is overplotted as a red line. Insert: Expanded 
plot of the distribution around the peak. Bottom: The observed distribution is compared to simulations derived for average magnetic fields 
of 100G, 200G and 300G respectively. 



